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Abstract 

We discuss minisuper space aspects a non empty Robertson- Walker uni- 
verse containing scalar matter field. The requirement that the Wheeler-DeWitt 
(WDW) operator be self adjoint is a key ingredient in constructing the physical 
Hilbert space and has non-trivial cosmological implications since it is related with 
the problem of time in quantum cosmology. Namely, if time is parametrized by 
matter fields we find two types of domains for the self adjoint WDW operator: 
a non trivial domain is comprised of zero current (Hartle-Hawking type) wave 
functions and is parametrized by two new parameters, whereas the domain of a 
self adjoint WDW operator acting on tunneling (Vilenkin type) wave functions 
is a single ray. On the other hand, if time is parametrized by the scale factor 
both types of wave functions give rise to non trivial domains for the self adjoint 
WDW operators, and no new parameters appear in them. 
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1 Introduction 



One of the simplest models of quantum cosmology is the Robert son- Walker (RW) 
minisuperspace. RW geometries describe homogeneous and isotropic universes. The 
RW geometry is defined by the line element 

ds 2 = -Nldrj 2 + a 2 ( V )dQ 2 3 . (1.1) 

In (1.1) the only dynamical degree of freedom is the scale factor a{rj). The lapse 
function N± is not dynamical , being a pure gauge variable. The quantity dVt\ is the 
standard line element on the unit three-sphere. We use units in which h — c—1 and 
G = M~ 2 = 3tt/4. 

The pure gravitational action corresponding to (|1.1|) is 

S « = TK~r I d i x^(R-2A) + -±- [ d 3 x^K 
167rG Jm 87rG JdM 

dnN±a 3 



a' 2 II 



a 2 \ A 



1.2 



N]_ J 3 

In (O) A is the cosmological constant^, M = I x S 3 is the space-time manifold, 
K = K\ is the trace of the second fundamental form of the space-like boundary 
dM = S 3 (defined by rj = const.) and the dot denotes differentiation with respect to 
rj. The Hamiltonian corresponding to ( |1.2|) is 

H = -N ± ^P 2 + a-g 2 a^ (1.3) 

where P a = dL/dd = —2ad/N± is the canonical momentum conjugate to 0(77) and 
g 2 = 4. It is assumed in ( (L3| ) that A > 0. Gauge invariance of ( |1.2j ) yields the 
Hamiltonian constraint 

-^-ri^-^^- (L4) 

The constraint ( |1.4|) requires a gauge fixing condition. A possible such gauge fixing 
condition is Nj_ = const. 7^ in which the time variable n becomes essentially the 

1 One can regard A as a pure cosmological constant or as the vacuum energy of some non trivial 
field configurations, or a sum of both. 
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proper time r. In this gauge the solution of the classical equations of motion (with 
initial conditions a(0) = g" 1 , d(0) = 0) is 

a ( T ) = S^cosh^r) , (1.5) 

which describes a universe that contracts from an infinite radius in the absolute past, 
reaches a minimum radius, a min = g~ x , and re-expands to infinity in the absolute 
future. 

Quantization of this simple system is accomplished straightforwardly in the coor- 
dinate representation by the usual operator realizations 

d 

a = a and P a = — i— . (1.6) 

da 

Neglecting operator ordering problems in the kinetic termQ the Hamiltonian con- 
straint becomes the Wheeler-DeWitt (WDW) equation |1|, for the wave function of 
the Universe: 

fid 2 \ 

+ a-g 2 a 3 *(o) = . (1.7) 



y 4a da 2 

Equation ( |1.7| ) is a Schrodinger equation = for a zero energy eigenstate of a 
mass m = 2 particle, moving in the one dimensional potential 



V(a) = a z -g z a\ (1.8) 

Despite the fact that the potential ( |1.8| ) is unbounded from below, it has the 
property that the time of flight of a classical particle from the largest turning point 
to infinity is finite, namely, 

dx 



<oo. (1.9) 



This may seem contradictory to ( |1.5| ) at first sight , because this equation implies 
r — > oo as a — > oo. However, this contradiction is only apparent, since what we call 
"time" in (|1.9|) is not the proper time r. In order that the classical particle that 



2 In this work we study only semiclassical solutions to the WDW equation, where operator ordering 
issues are not so important. 
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moves in the one dimensional potential ( |1.8| ) have the standard kinetic energy term 
(that is an implicit assumption in ( |1.9|) ), one must impose the "conformal time" gauge 
N± = a(t) on (|1.3|). Doing so, ( |1.3| ) becomes 

H = -^f + a 2 -g 2 a 4 y (1.10) 

The functional relation between the conformal time and proper time is 

f At 

t = J -— = tan" 1 ( smh(gr)) (1.11) 
and the classical solution ( |1.5| ) in the conformal gauge is 



a(i) = ^T (L12) 

It is clear from (|1.12|) that the "particle" reaches infinity indeed after finite conformal 
time, t = 7r/2. 

Eq. (|1.9|) suggests that one dimensional Schrodinger operators like (|1.10| ) are very 
similar to Schrodinger operators describing quantal systems defined on a finite seg- 
ment of the real line, despite the fact that they act on wave functions supported 
along the real positive half line. Quantal systems defined on a finite segment require 
boundary conditions as an essential part of their definition. In a similar manner, 
Schrodinger operators like (|1.10|) require boundary conditions on wave functions, en- 
suring probability conservation at infinity. Thus, such conditions extend ill defined 
Hamiltonians like (|1.10| ) into a self-adjoint form which generate unitary time evolution 
operators [|, 1 |, g 0. 

Considerations of extending the Hamiltonian ( |1.10| ) into a self-adjoint form seem 
irrelevant for quantum cosmological considerations at first sight. Indeed, solutions 
of the WDW equation ( |1.7| ) for the pure RW geometry ( |1.1| ) are always the zero 
energy eigenstate of (|1.10|) which spans a one dimensional Hilbert space on which 
( |1.10| ) is trivially self adjointQ This is, however, only an incorrect superficial state- 
ment. There are at least two important reasons to introduce self-adjoint extensions 

3 Note that the usual arguments that guarantee non-degeneracy of the discrete spectrum of one 
dimensional Schrodinger operators, and in particular-that the corresponding wave function are real 
up to an overall phase, are inapplicable for the operator in (1.1C) which is unbounded from below. 
We elaborate on this point in section 3. 
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of the WDW Hamiltonian of which ( |1.10[ ) is only a very simple case. First, we note 
that the null condition ( |1.4|) is just a special case of the general constraint enforc- 
ing reparametrization invariance on physical states in quantum gravity, namely, that 
they be annihilated by the WDW operator. This implies that one has to consider the 
WDW operator defined in a Hilbert space larger than the space of physical states, 
where one can apply it on various states and check whether they are physical or not. 
In this framework of constrained quantization it is always assumed that all relevant 
operators, and in particular, the constraint operators, are self-adjoint with respect to 
the Hilbert space inner product ||. In our particular case, this means that we have to 
include in the domain of definition of the operator ( |1.10| ) many eigenstates with non 



vanishing energy eigenvalues. In this case, the issue of a proper self-adjoint extension 
of (|1.10| ) becomes important. 



Second, note that when cases of a non-empty RW universes or perturbed RW 
universes |], [H| are considered, many non zero energy eigenstates of ( |1.10|) become 



relevant, even in the physical Hilbert subspace itself. This comes about because now 
the WDW constraint implies that the total Hamiltonian must annihilate physical 
states. To see this we observe that the total Hamiltonian may be written as 

Hwdw = H tot = Hfah — H (1-13) 



where — H Q is (|1.10|) and is the Hamiltonian of matter fields (denoted here by 



the field 0) and of gravitational perturbations (denoted by /^)Q- The minus sign in 
(|1.13|) results from the fact that the kinetic term of the conformal mode a in (|1.2|) 
appears with the "wrong" sign. In simple cases where a separation of variables applies, 
= \(pa)\X(f>,h) , the corresponding WDW equation reduces into the two equations 



and 



# O |0a> = E\(f> a ) (1.14) 



Hfrh\x<t>,h) = E \x<t>,h) (1-15) 



That is, the metric is g^ u = gffl + where is the RW metric in (1.1) 
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where E is some non-negative eigenvalue^] either of Hq or of H^h- Now, in principle, 
many eigenstates of H an d of (sharing the same eigenvalue) may enter into the 
quantum state of the Universe, especially if it turns out to be a quasi-classical state[|. 

Now that many different energy eigenstates of (|1.14|) and ( |1.15|) are used to con- 
struct physical states, one has to equip them with a time independent inner product 
in order to define the physical Hilbert space. We find that if we choose matter fields 
as our clock, self-adjointness of (|1.10|) turns out to be the necessary and sufficient 
condition for the existence of such an inner product. In this case there is a continu- 
ous two parameter family of non trivial domains of the self-adjoint WDW operator 
that are spanned by zero current (Hartle-Hawking type) wave functions, whereas the 
domain of a self adjoint WDW operator acting on tunneling (Vilenkin type) wave 
function is trivially a single ray. One can define a positive definite Hilbert space inner 
product using both kinds of wave functions. If, on the other hand, we choose the 
scale factor of the Universe to parametrize time, both wave functions give rise to 
non trivial domains, and are independent of any new parameters. However, positive 
norm physical states may be constructed using only Vilenkin type wave functions. 
Therefore, whenever we are free to choose either matter fields or the scale factor as a 
time coordinate, self-adjointness of the (spatial part of the) WDW operator dictates 
utterly different physical Hilbert spaces corresponding to different physical realities, 



and is therefore intimately related to the problem of time in quantum gravity|12 



In section 2 we consider a non-empty RW universe filled up with matter in the form 

of scalar fields where we derive ( |1.13| )-( [I.15| ) explicitly. We show how the requirement 

for self-adjointness of the (spatial part of the) WDW operator arises and point its 

relation to the problem of time. In section 3 we discuss self-adjoint extensions of one 

5 The requirement that E > is a consequence of the energy condition for matter fields. 
6 A quasi-classical state is a superposition of many energy eigenstates which are superimposed 
with some well behaved amplitudes, for example, 



^ qc (a) ~ J # B (a)exp 



E — E ^ 
AE 



dE 



where are solutions of (1-14 ) and Eq is the classical energy. 
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dimensional Schrodinger operators whose potential terms are unbounded from below, 
but satisfy (|1.9|) . In section 4 we apply such extensions to RW quantum cosmology 
and discuss the Hartle-Hawking and Vilenkin proposals for the wave function of the 
Universe in this context. Finally, we show in the appendix that the domain of the self 
adjoint WDW operator ( |1.10| ) acting on tunneling (Vilenkin type) wave functions is 
a single ray. 
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2 Non-Empty RW Minisuperspaces 



Consider a scalar field (j)(rj, x) which is coupled non-minimally to gravity. It is gov- 
erned by the action 

S m = - l -l d^^ld^d^ + ^ + V^)}-^ f <?xsf^)K<j)\ (2.1) 
2 Jm 1 J 2 JdM v 

where £ is the coupling constant. In the homogeneous and isotropic case, we have 
cf) = (j)(rj), and thus the total action (namely, the sum of ( |i.2p and (271)) is [TI[ 



Stot = S g + S m = I dr)N_ 
where % = 7ra 6 ^<p is the rescaled matter field and 



a 

N 2 



^3-12^2 

N 2 



+ U(a, X ) 



n<>.\) = « //V- ~ i " V(^ ? )-6£a 1 - 12 V 



The corresponding Hamiltonian is 

Hwdw = H to t = N± 

In the coordinate representation P a 
reads 



(2.2) 



(2.3) 



(2.4) 



— i-g- and P x 



the WDW equation 



1 d 2 



1 

+ T a 



12£-2 



d 2 



+ aU(a,x) 



(2.5) 



4 9a 2 4 «9x 2 

We concentrate on the simple conformally invariant case where V(4>) = and £ = 1/6, 
such that ( |2.4j) becomes 

' a U" 



WDW 



= H, 



tot 



\p^-\pi+" 2 -9 2 " i -x 2 ) 



(2.6) 



The WDW equation is therefore 
/ 1 <9 2 



V 49a 2 ' " * " ' 4<9 X 2 A ; " ■ (27) 

We solve ( |2.7| ) by separation of variables ^(a, x) = ipai^ipxix) i which results in the 
two coupled equations 

/ 1 d 2 x 



Ada 



2 + a 2 - g 2 a 4 ip a (a) 



1 d 2 



Ad X 2 



x 2 i> x (x) 



Eipa(a) 

EMx) ■ 



(2.8) 
(2.9) 



Quantization requires gauge fixing, namely a definition of "time". The freedom left 
in making such a gauge choice leads to the "problem of time" in quantum cosmology 
jn| |T3[. It is strongly related to the definition of a physical Hilbert space, and here 
is the place where considerations of self-adjointness come into playf|fl4| . One has to 
introduce an inner product in the space of solutions to the WDW equation ( |2.7| ) in 
order to define the physical Hilbert space. Defining the currents J^ 2 X = i^f* d a , x ^2, 



we have 



d V 1 - d J 



12 



0. 



(2.10) 



for any pair \1>2 of solutions of ( [2.7[ ) 

In our simple minisuperspace model we can choose the time variable either as the 
scale factor a or as the matter field \ W\ H3 - ^ we choose the scale factor a as 
our time, t = a, then 



<*l|*2>(a) =-i I d X 



^* 1 (a, X )da^2(a,x) 



(2.11) 



I a=t = const. 

is a natural inner product to use in constructing the physical Hilbert space. 

A well defined inner product among solutions of (|2.7|) is necessarily time indepen- 
dent in order that there be no conflict between time evolution of physical states and 
the definition of Hilbert space at each time slice t = constant. Deriving (|2.11 ) with 
respect to time we have 



(a) 



12 1 



-00 



(2.12) 



where we have used ( p.lOj ) and integrated over x- Therefore, time independence of 
(|2TTTp implies 

J^ 2 (-oo) = 0. (2.13) 



4 2 (+oo) 



This condition holds automatically due to finiteness of ( |2.11| ) which also leads to 
a discrete spectrum of (|2.9|) . The latter is the standard one dimensional harmonic 



oscillator spectrum, E n = u(n + l/2) = (n+1/2), for which J^ 2 (+oo) = Jl 2 (— 00) = 



7 The discussion in 1 14 is limited to the case of zero cosmological constant for which there is no 
need for self adjoint extensions of the WDW operator. 
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. The Schrodinger operator on the left hand side of ( |2.9| ) which is the spatial part of 
the WDW operator in (|2.6|) is clearly self-adjoint in this domain. Evidently, (|2.13| ) sets 
no restrictions on the Schrodinger operator on the left hand side of ( |2.8| ) and therefore 
self-adjointness of the WDW operator leads to no further consequences in this case. 
This implies, as we show in the next two sections, that a non trivial domain of the 
self-adjoint WDW operator includes either zero-current (zero-norm) wave functions 
or tunneling wave functions. 

Matters are utterly different when we choose x (and not a) as time. In this event, 
the roles of x an d a interchange. One should integrate ( [2.1 0|) over a, and the inner 
product is therefore 



(*i|* 2 ) (x) = -i da Wj(a, x) d x * 2 (a, x) 

•'0 L J \ x=t=con 

This inner product is time independent (x independent) only if 



(2.14) 



J^(+oo) - J^(0) = (2.15) 

for any pair of solutions of (|2.7|) . This condition is quite different from (|2.13|) . 
Unlike the latter, it is not satisfied automatically, because the potential energy in (|2.8|) 
is unbounded from below. Rather, it must be imposed, in order that the Schrodinger 
operator on the left hand side of ( |2.8j ) which is the spatial part of the WDW operator in 
( |2.6| ) for this choice of time parametrization, be well defined. Indeed, ( |2.15| ) is precisely 



the condition that the Hamiltonian on the left hand side of (|2.8| ) be symmetric with 
respect to this inner product. This can be accomplished by defining proper self- 
adjoint extensions of the Hamiltonian as we discuss in the next two sections. Namely, 
we show there that ( |2.15| ) implies that zero current wave functions (similar in form 
to the Hartle-Hawking wave function) give rise to a two parameter continuous family 
of non trivial domains for the self adjoint WDW operator whose spectrum {E n } is 
discrete and highly non linear in n. On the other hand, the domain of a self adjoint 
WDW operator acting on tunneling (Vilenkin type) wave functions is a single ray 
having a single eigenvalue. 
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Demanding that the Hamiltonian be self adjoint is equivalent to the requirement 
of unique time evolution of quantum states. We see that the two different choices of 
"time" in this model lead to different kinds of physical Hilbert spaces and therefore to 
two highly different physical "realities" . Thus, the requirement that the spatial part 
of the WDW operator be self-adjoint is intimately related with the problem of time 
in in quantum cosmology which manifests itself to its extreme in the model discussed 
here. 

We are not going to argue here which choice of time coordinate is superior to 
the other. However, we feel that the choice t = x has not received enough attention 
(at least as far as mathematical aspects are concerned), and we wish to fill this gap 
partially. On top of this, note that implications of self-adjointness of the (spatial part 
of the) WDW operator are less trivial here compared to those related with the other 
possible choice. 

In the non conformally invariant case, the WDW equation is generally non- 
separable and therefore much more complicated. However, in cases where V(x) has a 
narrow deep minimum at x = Xs we can assume to a first approximation that x does 
not fluctuate far away from that minimum and thus replace V(x) by its minimal value 
p s = V(xs)- In this approximation one thus merely shifts g 2 into g 2 ^ = g 2 + ir 2 p s . 

We close this section making some general remarks before turning to the formal 



discussion of the proper self-adjoint extensions of ( |2.8|) . Strictly speaking, the scale 
factor a in ( |1 . 1|) is defined on the ray a > 0. Therefore all Schrodinger operators 
involving a must be equipped with a suitable boundary condition at a = in order 
to make them self-adjoint. Such a boundary condition is a necessary datum purely 
from the Schrodinger theory point of view and for this reason we will impose such 
a condition on wave functions below. However, which boundary condition at a = 
must be imposed on (|1.10|) and (|2.8| ) in order to describe quantum cosmology is 
a highly controversial issue. Indeed, unlike its usefulness at large radii, the mini- 
superspace formalism we use in this paper might break completely at extremely small 
radii of the Universe because of the true singularity of ( |1 . 1| ) at the point a = 0. 
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Should this happen, ( |1.10| ) and (|2.8| ) will become useless at a —>■ as well as their 
solutions. The region a >> 1/g, on the other hand, is certainly in the validity domain 
of the minisuperspace quantization scheme. We may thus trust the wave functions 
resulting from the WDW equation only for large universe radii, as far as cosmological 
interpretations are concerned. 

As was discussed in the introduction, the fact that the time of flight (|1.9|) of a 
classical particle moving in the potential ( |1.8|) to infinity is finite has a very important 
consequence. Namely, it effectively turns the spectral problem involving the hamilto- 
nian on the left hand side of Ql.iO| ) into a problem defined on a finite segment, despite 
the fact that < a < oo. This calls for an appropriate boundary condition on wave 
functions at a = oo as well. The boundary condition imposed at a = must be con- 
sistent with the one set at a = oo. Thus, in principle, there is some influence by the 
a = endpoint (where the minisuperspace formalism is suspicious) on the asymptotic 
behavior of the wave function of the Universe as a — > oo (where the minisuperspace 
formalism is surely valid). We comment on this point in section 4. 
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3 Domain of the Self- Adjoint Hamiltonian Defined 

on the Ray < x < oo 



Consider a particle of mass m, moving in the one dimensional potential 

V(x) = x 2 — g 2 x A , x > 



(3.1) 



having energy < E < There are two classically allowed regions for this range 
of energies, x < x\ and x > x 2 . Here X2(E) > X\(E) > are the two classical turning 
points, namely, the two positive real roots of V(x) = E. 

A wave packet with energy distribution peaked at E will move in ( |3.1|) from 
X{ > X2 to infinity in a finite period of time 

dx 



t, 



(3.2) 



>2m[E - V(x)] 

raising the question what will happen to the wave packet as it "hits" the point at 
infinity or equivalently, how is probability conserved in such a system. Therefore, on 
account of (|3.2| ), unbounded motion in the potential ( |3.1| ) behaves in many respects 
as if it were bounded, and the point at infinity appears as if it were the end point of 
a finite segment [0. 

Probability conservation requires that the Hamiltonian governing this system be 
self-adjoint. This requirement on the domain of definition of the Hamiltonian is by no 
means trivial, since wave functions in the potential Q3.1| ) have only power-like decay 
while their first derivatives blow up at infinity, as can be most easily seen by writing 
down the leading WKB approximation 

1 



[E-V(x)]* 
+ C 2 (E) cos 



d(E) sin 



J s]2m{E - V(y))dy - j 



J pm{E - V(y))dy - j 



to a generic solution of the Schrodinger equation 

1 d 2 



2m dx 2 



+ V{x) 



V E {x) = E^ E (x) 



(3.3) 



(3.4) 
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in region x > x-i- Finiteness of in ( |3.2| ) means that ( |3.3| ) is square integrable 
for any value of E, but this is not true of its first derivative. Due to the fact that 
V(x — ► oo) — > — oo, local De-Broglie wave lengths of the particle become extremely 
short very quickly as it moves deeper into the classically allowed region rendering 
the WKB approximation more and more accurate as x — > oo Q. It is therefore 
enough to limit our discussion to the framework of the WKB approximation JF]. Such 
asymptotic behavior of and as x — > oo is in clear contrast with the 

exponential fall off of both bound state wave functions and their first derivatives 
in case of potentials that are bounded from below. In particular, it implies that 
there can be two square integrable linearly independent solutions of ( |3.4| ) sharing the 
same parameter E because their constant Wronskian need not vanish. Thus, square 
integrability is not sufficient to determine the spectrum, and what is needed is an 
explicit boundary condition at infinity which should be treated as if it were really a 
finite boundary point. 

For any two states and ^2 in the domain of the self-adjoint Hamiltonian we 
must have 

(i^il^a) = (*i|#*2>. (3.5) 



Using the coordinate representation and integrating by parts, (|3.5|) implies 

0. (3.6) 



1 ^2-*r 



doc doc f) 

This is precisely the consistency condition ( [2.15D encountered above. 

Considering non trivial domains of the self adjoint operator H, the current J 12 = 
\ {^2d x ^i — ^18x^2) must vanish at the two boundary points, since otherwise a 
particle reaching x = 00 will have to reappear at x = or vice versa. While such a 
periodic boundary condition is relevant for quantization of a particle in a finite rigid 
box|J, it is clearly improper here because in our case V(0) = while V(oo) = —00. 
We back this qualitative argument by an explicit calculation presented in the appendix 



8 Recall that this is also the region of scale factor a = x values where minisuperspace analysis is 
most valid anyway. 
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where we show that if the current J 12 did not vanish, the domain of the self adjoint 
operator H becomes trivial and collapses into a single ray. 

Therefore, studying non trivial domains, ( ^.6|) may be replaced by the stronger 
condition 



dx 



dx 



(x — ► 00 ) 



dx 1 dx 



(x = 0) = 0. 



(3.7) 



We show now that (3/7) will not hold (at x — > oo) for generic = and 
^2 = ^_b 2 ; -^i 7^ ^2, unless some special choice of the functions C\{E) and C 2 (E) 
in ( |3.3| ) is made. To make this point clearer, it is useful to introduce the phase 4> a (E) 
defined byM 



cot {<t> a {E)) 



C 2 (E) 
C^E) 



(3i 



where a is a parameter wave functions depend upon and will be determined later. In 
terms of 4> a (E), ( |3.3|) becomes 

Ci(E) 



[E- V{x)]±sm<f) a {E) 



-cos 



7T 



t 2m(E-V(y))dy + <p a (Ef) 

x 2 (E) 1 4 



(3.9) 



The cosine in ( |3.9| ) oscillates very rapidly as x — > oo and therefore ( |3.7|) will not be 
met (for Ei ^ E 2 ) unless the argument of the cosine becomes independent of E as 
x — > oo, namely, 



d_ 

dE 



x 2 (E) 



2m{E-V{y))dy + <t> a {E) 







x — > oo . 



(3.10) 



In order to solve ( |3.10| ) we need to specify an initial condition in E, this is how the 
parameter a gets in. One can choose 



t (E = a) = 



(3.11) 



and the solution of ( |3.10|) (subjected to ( |3.11|) ) is 



(E) = V2^\-r ^E-V{y)-s]a-V{y) dy + H - V(y)dy 

I Jx 2 {E) L J Jx9.(a) v 



rME) 

Ix2(a) 



(3.12) 



15 



where 22 (a) is the largest root of V(x) = a. Substituting the specific potential ( |3.1| ) 
into (|3.10|) and integrating over E we obtain 



\n+c 2 



; E < 

; E > 



(3.13) 



where z(E) = (1-AEg 2 )^, C a = £ {a) ( 2 lC{£=±)d( and /C(m) {£{m)) is a complete 
elliptic integral of the first (second) kind (we use S(m) below). Here we have also 
assumed a < so that C a < as well . 

From (|3~T^ ) we see that for x —* 00 ( |3.9|) approaches 



¥"\x 



00 



-cos 



X2(a) 



^2m(a-V(y))dy-^ 



(3.14) 



[£-F(:r)psin0 a (£) 

where the argument of the cosine is indeed E independent. It is clear now that ( |3.7| ) 
will hold at x —>■ 00 for any pair of functions tyj^ , of the form given by (|3.9|) 
and ( |3.12|) . These functions form a family of solutions of the Schrodinger equation 
(|3.4|) parametrized by the single (real) parameter a. These are not eigenstates of the 
Hamiltonian, as we have not imposed (|3.7| ) at x = yet. To this end we note that a 
necessary and a sufficient condition for vanishing of the probability current at x = 
is 

*' E (x = 0) 



*e(x = 0) 13 
where f3 is a fixed real number || |], ||, 

In a similar manner to ( |3.9| ) we write the WKB solution of ( |3.4j ) as 



(3.15) 



¥£\x < Xi) 



cm 



-cos 



[E -V(x)]*sm£p(E) 

and ( |3.15| ) fixes 

£,p(E) = arctan 



xi(E) 



2m(E-V(y))dy + t l3 (E) 



xi(E) 



yj2m(E-V(y))dy. 



(3.16) 



(3.17) 
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Eigenstates of (p.4|) are obtained by matching ( |3.9| ) (subjected to (|3.12| )) and 
(|3.16|) - (|3.17|) . Using the ordinary WKB matching conditions we obtain the corre- 
sponding "Bohr-Sommerfeld" quantization condition on E, namely, 



where 



A(E) 



4tan0 a (£) = e" 2A{jB) tan 



x 2 (E) 



xi(E) 



2m\E - V(y)\dy 



(3.18) 



/ 2mx 2 (E) 



8 1 



xi(Ef 
x 2 (E)< 



-2g'x((E))C[l- 



Xi{Ef 
x 2 (E)< 



. (3.19) 



So far we have concentrated on energy range < E < Our discussion may be 
extended in a straightforward manner to the complementary energy ranges E > 
and E < as well. For example, the quantization condition corresponding to E > 
reads 



tan 



^2m(E - V(y)) - ^2m(a - V(y))) dy 







(3.20) 



V2mE 

Note the explicit dependence of energy eigenvalues and eigenf unctions upon a and 
(3. The spectrum is therefore a two parameter family {^^'^(x)} , parametrized by a 
and j3 as it should be in this case of separeted boundary conditions according to the 
general theory of self-adjoint extensions^, |], ||, |?J. Due to Q3.2| ) the point at infinity 
appears as if it were a finite endpoint and the whole quantum system behaves as if 
it were defined on a finite segment, where a and (3 parametrize boundary conditions 
at the two endpoints. The set of functions (x)} spans the domain of the self- 

adjoint Hamiltonian, namely, the space of all square integrable functions that satisfy 



Note that the WKB density of states 



, (E) — — [°° dx 

t Jo J2m[E - V{x)] 



(3.21) 
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which is proportional to ( |3.2| ) is finite (as long as E ^ ^r)f\- The spectrum must be 
therefore discrete for any E ^ (with an accumulation point at ^). It is moreover 
bounded neither from below, nor from above and contains therefore an infinite amount 
of discrete states. Thus, the domain of the self-adjoint Hamiltonian is the set of 
all discrete linear combinations of the form 




(3.22) 



n 



where J2 n \ c n\ 2 < oo which yields an infinite dimensional Hilbert space. 



For E — ► ( 3.21 ) diverges logarithmically in \E — 5=3- 1. 
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4 Quantum Cosmological Implications 

The most important observation made in section 2 (as far as the simple cosmological 
model discussed there is concerned) is that time independence of the inner product in 
the definition of the physical Hilbert space leads to the requirement that the spatial 
part of the WDW operator, which is a one dimensional Schrodinger operator, be self 
adjoint. We saw in that section that one can choose time parametrization either in 
terms of the matter field \ or in terms of the scale factor a of the Universe and that 
these two choices of time parametrization lead to different physical realities. In our 
view this is an extreme manifestation of the problem of time in quantum cosmology. 

Our discussion in the previous two sections makes it clear that one may trace this 
discrepancy of physical realities to the fact that the two choices of time parametriza- 
tion lead to two utterly different physical Hilbert spaces. The reason for this difference 
stems directly from the requirement that the (spatial part of the) WDW operator be 
self-adjoint. In this section we sharpen this distinction and investigate its cosmologi- 
cal implications in more detail. 

Let us first concentrate on parametrization of time in terms of the scale factor a. 
From the discussion following (|2.11 ) it is clear that a generic physical state has the 



form 

X) = J2 C n^an(a)lp xn (x) (4-1) 
n>0 

where c n are complex constants, ip xn {x) is the normalized harmonic oscillator eigen- 
state corresponding to eigenvalue E n = (n + |) and ip an (a) is a normalized solution 
of ( |2.8| ) with parameter E = E n so that matter excitations are seemingly those of 
a free field. Because of orthonormality of the ip xn {x) the inner product ([2.11|) of any 
two such states is simply 

^ (1) l^ (2) )(.) = E^i 2r ci 1) (4.2) 

n>0 



10 Recall from our discussion in section 3 that these solutions are normalizable because of the 
finiteness of (3.2), as long as E ^ -^p. 
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where j n = —i[^)* an {a)ip an {a)' — ip* an {a)' 'Vw(a)] is the "Schrodinger" current carried 
by ip an (a). This current is clearly a independent, making (|4.2|) time independent in 
accordance with ( [2.12| ). The "Schrodinger" current carried by ip a n{ a ) must be positive 
in order that ( f4.2| ) be positive definite. This positivity condition becomes simply 



|5i(£)| 2 - \B 2 (E)\ 2 >0 



(4.3) 



where B\ and B 2 are the amplitudes of the outgoing and of the incoming waves in 
the generic solution of ( f2~8[ ) to the right of the largest turning point a 2 (V(a 2 ) = 0)0 

1 



ipaE{a > a 2 ) 



[E - V(a)} 
B 2 (E)exp 



5i(E)exp 



a 2 (E) 



yj4(E-V(y))dy 



+ 



il ^A(E-V(y))dy 

o 2 (b) 



(4.4) 



For B 2 = 0, ( |4.4j ) gives pure expansion modes of the metric. The lowest mode 
(corresponding to n — 0) is the original wave function suggested by Vilenkin ||1 8|| . 
Wave functions of the form (fO|) subjected to ( |4.3|) may be referred to as "generalized 
tunneling wave functions" [|nj. Note further that if we have \Bi(E)\ = \B 2 (E)\ for 
all ipan{a) in (EO) the latter describes standing wave modes of the form 

B(E) 



^ aE {a > a 2 ) 



cos 



7T 



4(tf-V(y))dy + 0(£)-7 
02 (JS) 4 



(4.5) 



[E-V{a)\* I 

carrying no current, that is j n = for all n. The lowest of these modes (n = 0, <f)(E) = 

0) is the wave function suggested by Hartle and Hawking ||T9|, p0| , |TT| . Following (|4.2|) 

one cannot associate probabilistic interpretation with such wave functions because all 

physical states in fl4.1| ) become zero norm states. 

We turn now to the case where time is parametrized by the matter field x — vra0. 

Here the spatial part of the WDW operator is the Schrodinger operator in ( |2.8|) . 

which requires non-trivial self-adjoint extensions as we saw in the previous section. 

It is shown in the appendix that in this case the domain of the self adjoint WDW 

11 Recall that here E — E n > 0. For all practical uses we assume further that E < so that 
a 2 > 0. 
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operator acting on tunneling (Vilenkin type) wave functions is trivial, namely, a single 
ray and the spectrum shrinks to a single point E = E . It is very interesting that 
the seemingly innocent requirement that the WDW operator be self adjoint is so 
powerful that it allows only a single physical state to exist. We stress that such a 
trivial domain occurs only for tunneling type wave functions. This may have the 
far reaching cosmological implication as the mechanism that selects the unique wave 
function of the Universe[f[ 

The zero current wave functions span non-trivial domains of this self adjoint WDW 
operator. In the classically allowed region a > ci2 (which is the region where the 
minisuperspace approach is most reliable in the first place) the WKB approximation 



to the eigenf unctions spanning the non trivial domain are found from (|3.9|) as 



^\a>a 2 )= ° {En) (cos 



f JA(E-V(y))dy + MEn)-l 

Ja 2 (E n ) V 4 



(4.6) 

where 4> a {E) is given by ( |3.12|) and E n = E^'$> are solutions of ( |3.18|) . The modes 



given by (|4.6| ) that appear in physical cosmological states must have non-negative 
energies E n > otherwise the energy condition for matter fields will be violated. We 
note here that for most practical purposes only energies in the range < E n < 
are important. 

Clearly there is a continuous family of such domains parametrized by the two real 
continuous variables a and (3. The functions ( |4.6| ) are standing wave modes of the 
geometry which are similar in form to the Hartle-Hawking wave function, but do not 
correspond to the "no-boundary" proposal for E > 0. 

A generic physical state (corresponding to time parametrization by the matter 
field x = has the form 

* M Kx) = E^ ) («)^W (4-7) 

n>0 

where now ip xn {x) are solutions of ( |2.9|) with parameter E = E n . The functions 

12 Note that the requirement of self adjointness by it self cannot fix E which will hopefully be 
determined by the exact short distance Hamiltonian of quantum gravity. 
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ip X n{x) are therefore generic parabolic cylinder functions |T7|, |21[| rather than harmonic 
oscillator eigenstates because the {E n } in (|3.18| ) are obviously not harmonic oscillator 
eigenvalues. It is convenient to express the general solution of the parabolic cylinder 
equation (|2.9|) in terms of the two linearly independent functions U n and V n asfl~7| 



^xn(x) = hU n (2 X ) + c n V n (2 X ) , (4.8) 



where 



U n (z) = D E i(z) and 



VJz) 



T(h - E n 



• 2 



D E i(-z) - sm(itE n )D E i(z) 



TT 

Here D\(z) are Whittaker functions, and {b n , c n } are complex constants. 

The asymptotic behavior of the particular combinations of Whittaker functions 
appearing in ( |4.9|) as z — > ±00 is[fT^, pi] ?7(z) ~ e^^ 2 1^;^^ - i and V(z) ~ e ±2;2 |2| =F ' E ~^, 
and their Wronskian is \ -. 

V 7T 

Because of orthonormality of the V'an ( a ) the (x independent) inner product 
Q2.14 ) of any two states of the form (|4.7|) becomes 



(*U\^\ X) = -if- - c^Hf) . (4.9) 

V 71 n>0 

In order that physical states have positive norm we must impose Im(b* n c n ) > 0, a 
possible choice being c n = z6 n Q- 

The physical Hilbert space inner product is ( fOf) and it is finite provided {b n , c n } 
are properly restricted, for example, by demanding that J2n(\b n \ 2 + l c n| 2 ) < 00 ■ Fol- 
lowing (|3.22|) , note however, that in order that ( |4.7|) as a whole be in the Hilbert 
space of the self adjoint operator ( ^.8| ) we must have N(x) — J2 n l^x^ix)^ < 00 • 
This condition will hold generally for finite values of x provided ( |4.9| ) is finite. How- 
ever, from the asymptotic behavior of U n and V n we know that N(x) blows up like 
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Note that using harmonic oscillator eigenstates for the ^ xn {x) makes a ll s tates in (4.7) zero 



norm states in the same way that using (h£) for the V'an(a) renders all states (4J) zero norm states 
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e x as |x| — > oo throwing ( (4.7|) out of the domain of (|2.8|) at |x| = oo, event hough 
(|4.9|) remains finite. The limit \x\ — > oo is easily attainable for any finite value of 
the scalar field as the scale factor a blows up to infinity because x = na<ft- This 
apparent difficulty deserves further investigation, but we are not going to do so in 
this paper. We only mention that this problem can be avoided in the single tunneling 
wave function domain, because the single energy Eq can be always chosen as one of 
the harmonic oscillator eigenvalues with the corresponding square integrable ip X Eo{x)- 

We see that a single (and relatively simple) physical system, namely, a conformal 
scalar field coupled to gravity, is described by two different and incompatible physical 
Hilbert spaces, corresponding to the two different possible time parametrizations. 
This implies that two highly different physical "realities" correspond to the same 
field theory coupled to gravity which in our opinion is an extreme manifestation of 
the problem of time in quantum gravity^. 

In the model discussed above, time parametrization by the matter field \ forces us 
to use the gravitational wave functions fp^'^{a) in ( f4.6| ) (for a > g^) in constructing 
physical states ( |4.7|) in the non trivial domain. This means that there is a continuum 
of distinct physical Hilbert spaces that are parametrized by the two real variables a 
and (3 which are the domains for the non trivial self adjoint extensions of fl2.8|) . Thus, 
specifying a physical state requires fixing a and (3 first. 

The most urging physical question that arises concerns the cosmological interpre- 
tation of a and (3: is there really a continuum of distinct WDW operators or are these 
parameters fixed somehow by a yet unknown dynamical mechanism associated with 
quantum gravity^]? If there is such a continuum of WDW operators, which of them 
corresponds to the "real" Universe? 

To partially answer these questions, recall that the wave function of the Universe 
is a unique solution of the WDW equation that must be singled out of all the other 
mathematically possible solutions whose general WKB form is fl4.7|) . One then gener- 



ally argues [13] that the distinguished semi classical solution which is the wave function 



14 Using the terminology of Kuchaf fil|, this is the "Multiple Choice Problem". 
15 As happens, for example, with the 9 angle in the Standard Model. 
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of the Universe will be picked up by matching it to a yet unknown solution of the 
WDW equation corresponding to the exact theory of quantum gravity that governs 
small geometries. In the case discussed here more is required, namely, that such a 
matching condition will teach us something about a and (3 as well. It might be that 
the parameter (3 is an artifact of our minisuperspace approach, because it defines the 
boundary condition ( |3.15|) at a = where the RW minisuperspace approach as a 
whole probably breaks down. Thus, it is plausible that if we knew the exact WDW 
operator for small geometries the ambiguity associated with (3 would be lifted either 
by showing that (3 is indeed a minisuperspace artifact, or by fixing its value if it has 
anything to do with the exact Hamiltonian. Unfortunately, we cannot see how such 
considerations apply as far as the parameter a is concerned. Indeed, a specifies the 
boundary condition at a — > oo where the semiclassical approximation is perfectly 
valid, and moreover, it is completely independent of /3. It seems to us highly unlikely 
that small geometry quantum gravity effects will be able to remove the ambiguity in 
the WDW operator and the wave function of the Universe associated with a. 

As a matter of fact, a has the following interesting cosmological interpretation: 
The original Hartle-Hawking wave function of the Universe [H], |20|, |TT] \1/ hh{o- > 0,2) ~ 
(|U(a)|) -1 / 4 cos (ja 2 (E) \l^\V(y)\dy — is obtained from the semiclassical approx- 
imation to the Euclidean path integral based on the "no boundary" proposal p2f. It 
is governed by one of the two Euclidean solutions to the classical equations of motion 



having the larger action. It has been shown, however, in [p3] that the contribution to 
the path integral from the other classical solution is relevant as well. It follows than, 
that if one considers contributions to the path integral weighed by real coefficients, 
one obtains a generalization of the Hartle-Hawking wave function given by ( f4.6| ) eval- 
uated at E — 0. Thus, cos0 a (O) weighs the contribution of the larger action saddle 
point and sin0 a (O) weighs the other one. 

The "no boundary" proposal of Hartle and Hawking points at the special role 
played by zero energy E = 0. It corresponds to regular Euclidean geometries, whereas 
singular Lorentzian geometries occur for E > 0. Note further that (3, whose inter- 
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pretation is problematic from cosmological point of view, disappears from ( |3. 17] ) and 



(|3.18|) at E = 0. It is therefore interesting to study the consequences of including 



E = in the discrete spectrum of (|3~TB1 ). Using ([TTJ), Q3T7D and (^T|) at £ = 



and demanding that ( p,18|) be satisfied at E = as well, we obtain the following 
functional relation between g 2 and a 

where C a was defined following ( |3.13| ). Thus, for a given value of a ( |4.10| ) quantizes 
the cosmological constant A = 3g 2 . Note further that Q4.10|) has a maximal positive 



solution (yf^ aa , which sets an upper bound for the cosmological constant. For realistic 
values of the cosmological constant g 2 « 1 Eq. (|4.10| ) reduces to 



A n =^i , n»l^l. (4.11) 
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We close this section by making some comments regarding possible extension of 
our work. For simplicity, we have focused our discussion on RW cosmologies that are 
isotropic and homogeneous. Generalization of our work to more complicated cosmolo- 
gies (presumably with positive cosmological constant) is clearly required. Moreover, 
implications of uniqueness of the tunneling wave function as the wave function of 
the Universe in case of time parametrization by matter fields to inflationary models 
should be clarified. Work on these points is in progress. 
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A Appendix 

Consider the self adjoint Hamiltonian H = ^P 2 +V(x) ; x > and its domain V. We 
assume that V(x) is unbounded from below such that (1.9) holds. For simplicity we 
also assume that V(0) = is a local extremum of V(a;)0. We show in this appendix 
(within the framework of the WKB approximation) that if T> consists purely of right 
moving (or of left moving) waves, it is trivial, namely, a single ray. We also present 
indications that this is the case for generalized tunneling wave functions as well. 

We first consider right moving waves, but the results are the same for left moving 
waves. Let x > (E) be the largest turning point of a classical particle with energy E 
moving in a general potential V(x) satisfying (1.9). A right moving WKB solution 
to the right of x > (E), is 

**"' (X > *>> = [2m(E -V( x) )f^ (' jCw VMB " m)dy) . (A.l) 
where A E is a complex amplitude. 

For energies above the maximum V max of V(x) (that is E > m the case of 
(3.1)) there are no turning points and ( [A.l]) is the wave function for all x > so that 
x > = 0. 

Consider the current 

J 12 (x) = l -{* 2 d x n - nW(x) . (A.2) 



At x = we have 



J 12 (0)) =\a*a 



F~>V 2' B l 2 

Ej j> vmax i-i 



eJ \e x . 



(A.3) 



while for x — > oo the current becomes 

J 12 {x -> oo) — ► A* Ei A E2 e- i5l2{oo) . (A.4) 

Here 6 12 (x) = J^ >(a) pm{E 2 - V(y))dy - f* >{El) pm{E x -V{y))dy. From (Q 
and ( |A.4| ) we see that J 12 (0) = J 12 (oo) if and only if E\ = E 2 , so T> contains at most 
only a single eigenstate in the energy range E > V max . 



Hi 



Note further that V"(0) > is a necessary condition for E ~ to be in the spectrum. 
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For energies below the lowest local minimum V m i n of V(x) (that is, E < in the 
case of (3.1)) there is only one turning point, x (E), and the current at x = is 

1 



(J u (x = 0)) 1 



2 -A* A E% |cosh(a - 6) 

_2 fl e -(?i+6) _ e 6+£2 
U 



^\ 1/4 /^\ 1/4 ' 

£ 2 \ 1/4 /^i\ 1/4 



(A.5) 



Si 



£/2 



where = Jq ^ \/ 2m(V (y) — Ei)dy. The norm squared of the current at x = is 
I J 12 ' 



2 



A* cosh^-6) 2 




+ 



_ e -(6+6) _ e 5i+6 




+ 



Ei 
E2 

E 2 





while from (|A.4j ) the norm squared at x — > 00 is 



J 12 (x 



00 



\A* A, 



(A.6) 



(A.7) 



The last two expressions are equal if and only if E\ = E2 and thus also in this energy 
range T> contains at most a single eigenstate. 

For energies in the range V m in < E < V max there are more turning points and 
J 12 (0) is sensitive to the details of the potential, but one can show that the domain 
T> is trivial as well. For concreteness let us consider the potential (3.1). In the energy 
range < E < there are two (positive) turning points, < x < < x > . Using the 
WKB matching conditions, the wave function to the left of x K (E) which corresponds 
to (TA~T| ) is 

( e -A(E) 



(WKB) 



[X < X, 



-A E e in / 4 



[2m{E-V{x))Y/* 



-sin[0 E (x)] + 2ie A(E) cos[(p E {x)} 



(A.f 



where <j> B (x) = J*< (E) ^2m(E - V(y))dy - f and A(E) is given by ( gig ). The 
current at x = becomes 

J 12 (0) = ^A* E ^A E2 I [(ci2 + c 12 1 )cos0icos0 2 
+ (di2 + c?i 2 1 )sin0isin0 2 

— i (a i2 + a^ 2 1 )cos0isin0 2 - (&12 + fr^sin^icos^ } , (A.9) 
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where 



«12 



Cl2 



-(A(E 1 )+A(E 2 )) 

4 V£ 2 7 



£i\ 1/4 



E-> 



-(A(Ei)-A(B 2 )) 



&12 = 
C?12 



1 

4 V^i 

' E 2 



1/4 



1/4 



-(A(E 1 )+A(B 2 



}) (A.10) 



-(A(Ei)-A(£ 2 )) 



and 



(0) 



/ 2mx 2 (E i 



5 



1 - 4£^ 2 /C 



7 • A.11 



The norm squared of ( A.9|) can be written in the form 

1 

4 

+ 



I J 



12 



-U:. .4. I 2 {(ci2 + c^fcos 2 ^ + (d la + d^sra^i 



E 2 



a-12 + ai 2 1 )cos0isin0 2 - (&12 + 6i 2 1 )sin0icos^2 
C12 + c5" 2 1 )cos0isin0 2 - (^12 + ti^ 2 1 )sin0 1 cos02 



(A.12) 



For fixed E\ one may consider the right hand side of ( A.12|) as a function of E 2 . Using 
( |A.ll| ),( fA.ll| ) and ( |3.19| ) it can be shown that this function has a global minimum at 
E 2 = Ei where trivially J 11 (0) = J 11 (00). Thus, for any other E 2 7^ E\ (in this 
energy range) we have \J 12 ' 



> |J 12 (oo)| and (3J5) does not hold. Therefore T> 
contains a single eigenstate in this E range as well. 

So far we saw that T> contains at most a single eigenstate in each of the energy 
ranges E < V m i n , Vmin < E < V max and E > V max . Considering now self adjointness 
of the Schrodinger operator in this finite dimensional domain T> we easily see that T> 
is actually only one dimensional (because otherwise | J 12 (0)| > | J 12 (oo)| for any pair 
of different such states) and the spectrum of H shrinks to a single point E = Eq. 

For generalized tunneling wave functions which are not purely right (or left) mov- 
ing waves, the explicit form of the (non-zero) J 12 (0) is more complicated and we 
currently do not have a complete proof that T> is one dimensional, but we present 
below indications that it is finite dimensional at most. The WKB wave function to 
the right of the largest turning point is 

1 



(WKB) 

V E (x>x > , 



(e - v(x)y/ 4 



(e) 



2m(E - V(y))dy 
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+ B E exp (-ij* e yj2m(E-V(y))dy 



(A.13) 

l X> (E) 

where A E and B E are complex amplitudes. For self adjoint Hamiltonians with po- 
tentials satisfying (1.9) that are defined over an infinite dimensional domain T> the 
spectrum is necessarily unbounded from below and from above. In this case one 
can choose two eigenstates in V with energies E x and E 2 above V max , such that 
E 2 /Ei — ► oo. To leading order in (Ei/E 2 )^ « 1 the current at x = is 

1 / E \ 

J 12 (0) ~ - ( — ] ( A* A„ — B* B„ + B* A„ — A* B„ ") . (A.14) 
On the other hand, the current at infinity is 



ji2 (a . ^ oo) ^ Al A Eo e~ i5 ™ - Bl B E2 e^. (A. 15) 



A generalized tunneling wave function is a wave function for which \A E \ > \B E \. For 
such wave functions the norm of (|A.14 ) diverges as E 2 /Ei — > oo, while the norm 



of ([A. 15 ) is finite. Thus, J 12 (0) ^ J 12 (oo) and therefore, contrary to our initial 



assumption, the spectrum has to be bounded from above. Similarly, using ( A.6 ) for 
negative energies we find that J 12 (0) diverges when |J5 2 |/|J5i| — > oo, and the spectrum 
must be bounded from below as well. There must be a finite minimum gap between 
energy levels in the band < E < and thus the number of eigenstates is finite 
at most. Because this argument involves the limit E 2 — > oo it is insensitive to the 
details of V(x) near x = as long as V(x) and V'(x) are regular at the origin. 

J 12 (x = 0) can be equal to J 12 (x — > oo) for E 2 jE\ » 1 only if they are both 
zero. Only in this case J 12 (x = 0) will not diverge as E 2 /E\ — > oo, and we can have a 
non-trivial {infinite dimensional) domain for the self-adjoint Hamiltonian. This was 
studied in section 3. 
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